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ABSTRACT 


The  degradation  in  the  sensitivity  of  a FLIR  as  a result  of 
the  absorptance  and  emittance  of  the  infrared  window  material  is 
considered.  An  analytical  treatment  of  the  problem  is  presented 
along  with  some  calculated  results  for  the  PAVE  TACK  system  and 
three  possible  window  materials. 
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I - INTRODUCTION 


The  PAVE  TACK  system  currently  utilizes  a zinc  sulfide  window 
in  conjunction  with  the  AAQ-9  FLIR.  There  is  some  concern  and 
speculation  as  to  how  much  loss  in  sensitivity  of  the  system  is 
attributable  to  the  window,  and  how  much  the  sensitivity  would  be 
gained  by  the  use  of  an  improved  window  material  such  as  gallium 
arsenide  or  zinc  selenide.  This  report  establishes  a method  of 
isolating  the  loss  in  system  sensitivity  due  to  the  window  and 
presents  results  for  several  window  materials  and  conditions. 

The  theory  behind  this  report  is  not  new  and  calculations  are 
similar  to  those  made  before.  The  uniqueness  of  the  report  is  that 
effects  of  atmospheric  transmittance/radiance  are  included  as  well 
as  specific  and  current  data  on  AAQ-9  system  response  and  properties 
of  three-window  materials.  The  established  methodology  simplifies 
new  calculations  whenever  new  data  becomes  available  or  when 
additional  trade  studies  are  appropriate. 

It  should  be  noted  that  there  are  many  factors  which  go  into 
the  selection  of  the  best  window  material  for  a particular  applica- 
tion. The  loss  in  sensitivity  due  to  a window  must  be  considered 
along  with  strength,  hardness,  erosion  resistance,  solubility, 
coating  adhesion,  homogeneity,  cost,  producibility , etc.  in  making 
this  selection. 

This  report  was  prepared  by  Douglas  W.  Amlin  with  considerable 
assistance  from  Ronald  T.  Vantrease  and  Roberto  F.  Soto,  all  personnel 
of  the  Imaging  Systems  Branch,  Mission  Avionics  Division,  Directorate 
of  Avionics  Engineering,  Deputy  for  Engineering,  Wright-Patterson 
Air  Force  Base,  Ohio. 
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II  - PURPOSE 


The  purpose  of  this  paper  is  to  calculate  Sensitivity 
Reduction  Factors  (SRFs)  for  three  different  window  materials; 

ZnS,  ZnSe  and  GaAs.  The  PAVE  TACK  FLIR  system  with  known  relative 
spectral  response  in  the  7.5  - 12.0  micron  region  is  assumed. 
Results  are  presented  which  show  how  the  SRF  for  each  window 
material  varies  as  a function  of  window  temperature.  In  addition, 
atmospheric  effects  on  SRF  are  presented  based  on  atmospheric 
transmission  and  radiance  factors  computed  for  a variety  of 
conditions. 


Most  airborne  FLIR  systems  are  used  in  conjunction  with  an 
infrared  window  which  acts  as  a barrier  between  the  FLIR  and  the 
dynamic  environment  in  which  the  aircraft  is  operating.  The 
infrared  window  is  generally  not  perfectly  transparent  throughout 
the  spectral  bandpass  of  the  FLIR  and,  consequently,  affects  the 
performance  of  the  FLIR  by  a combination  of  absorption  and 
emission  of  detectable  radiation.  It  can  be  shown  that  the  loss 
in  thermal  sensitivity  induced  by  an  IR  window  is  proportional  to  the 
ratio  of  the  Signal/Noise  Ratio  (SNR)  with  and  without  the  window, 
(SNR)W  / ( SNR ) 0 or  Sensitivity  Reduction  Factor  (SRF).  After 
making  a few  general  assumptions,  this  factor  can  be  calculated 
without  further  consideration  of  specific  features  of  the  sensor 
system  under  consideration. 

It  has  been  shown  by  Klein1*2  that  a "simple  and  elegant 
formalism  exists  for  assessing  the  degradation  in  signal  to  noise 
ratio  resulting  from  the  presence  of  a partially  transparent 
window  at  the  entrance  aperture  of  a FLIR  sensor."  The  degrada- 
tion involves  a reduction  in  target  signal  and  an  increase  in 
system  noise  or: 


HEFFW/HEFFC 

NEIw/NEIq 


(1) 


HEFFW/HEFF0  is  the  ratio  of  effective  signal  irradiance  with  and 
without  the  window  and  is  calculated  from: 


heffv  _ W»(Ts)  a>. 

HEFFq  " Wx  (Tg)  <n 


NEIW/NEI0  is  the  ratio  of  noise  equivalent  irradiance  with  and 
without  the  window  and  is  calculated  from: 
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The  individual  terms  in  these  equations  are  defined  as  follows: 
T(X,TW) 


eU.t„) 


MV 

MV 


Transmittance  of  the  window  at  a particular 
wavelength  and  temperature. 

Blackbody  spectral  radiant  exitance  evaluated 
at  the  temperature  of  the  background. 

Emittance  of  the  window  at  a particular  wave- 
length and  temperature. 

Spectral  quantum  exitance  function  evaluated  at 
the  temperature  of  the  background. 

Spectral  quantum  exitance  function  evaluated 
at  the  temperature  of  the  window. 


The  quantity  SNI^/SNRq  may  properly  be  defined  here  as  a 
sensitivity  reduction  factor  (SRF)  since  the  noise  equivalent 
temperature  (NET)  of  a system  with  a window  can  be  defined  as: 


or 


SNRq 

NETv  = SNR^  NETo 


NETW  = 1/SRF  NET0 


U) 


(5) 


Likewise,  the  minimum  resolvable  temperature  (MRT)  of  a 
system  at  target  spatial  frequency  (f^,)  with  a window  can  be 
calculated  from: 


MRT^f T ) = ?w(fT)  MRTQ(fT) 


(6) 


SRF 


where  represents  the  window-related  contribution  to  the  modula- 
tion transfer  function  of  the  system  and  is  obtained  from: 


^ ,r  s MTFw(fT) 


MTFQ(fT) 


(7) 


Based  on  these  equations , it  is  seen  that  the  effect  of  a 
window  on  the  sensitivity  of  a FLIR  system  can  be  reduced  to  a 


U 
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single  figure  of  merit: 

SRF  = SNR  /SNR  (8) 

w o 

This  value  can  be  calculated  with  very  little  Knowledge  of  the 
specific  features  of  the  system.  However,  it  should  be  noted 
that  the  foregoing  assumes  (l)  that  the  window  is  at  a uniform 
temperature  and  exhibits  lambertian  characteristics,  (2)  the 
detector  package  includes  a cold  shield  which  blocks  all  internal 
noise,  (3)  no  absorption  by  the  sensor  optics  or  atmosphere  and 
(k)  a flat  detector  response  throughout  the  spectral  bandpass  of 
the  system. 
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IV  - FORMULATION 


In  order  to  utilize  the  basic  theory  and  equations  of  Klien 
to  evaluate  window  degradation  for  a practical  situation,  two 
modifications  can  be  made.  First,  a detector  factor  can  be  worked 
into  the  equations  and  secondly,  the  effect  of  the  atmosphere  can 
be  taken  into  account.  Without  proper  inclusion  of  these  consid- 
erations, it  is  impossible  to  precisely  isolate  the  effect  of  the 
window  on  system  performance. 

Inclusion  of  a detector  factor  into  the  equations  presents 
little  problem  if  the  system  response  is  known.  A detector  factor, 
D(\),  must  be  inserted  as  a factor  in  each  term  of  the  equations 
for  HEFFW/HEFF0  and  NEIW/NEI0,  and  must  include  the  relative 
response  of  the  detector,  the  transmission  of  the  system  optics, 
and  the  transmission  of  the  bandpass  filter.  Setting  all  values 
of  D( X ) equal  to  1.0  is  equivalent  to  disregarding  the  factor  or 
"backing  it  out"  of  the  calculation.  Setting  peripheral  values 
of  D( X ) equal  to  zero  is  equivalent  to  setting  the  bandpass  of  the 
system. 

Inclusion  of  atmospheric  effects  into  the  calculation  is  a 
bit  more  complicated.  The  effect  of  the  atmosphere  is  analogous 
to  the  effect  of  the  window;  i.e.,  a loss  of  signal  and  an  increase 
in  noise.  Values  for  atmospheric  transmission,  A(x),  and  values 
of  atmospheric  radiance,  R^l1^),  can  Le  calculated  by  use  of  the 
LOWTRAN  4 computer  program.3  In  this  program  Rx(Tb)  includes  the 
emission  from  the  boundary  (earth)  at  a specified  temperature  plus 
emission  of  the  atmosphere  for  the  specified  slant  range  and 
weather  conditions.  Since  the  boundary  is  included  in  Rx(Tn),  the 
quantity  Qx(Tb)  in  Equation  (3)  is  no  longer  needed  when  RX(Tb)  is 
used.  However,  the  units  of  RX(TB)  must  be  changed  from  WATTS/cra2- 
STER-micron  to  PHOTONS/cm2-SEC-micron,  RQx(Tfi),  to  allow  combina- 
tion with  Q\(TW),  the  noise  coming  from  the  window.  Recalling 
that  QX(T)  = WX(T)  X/hc,  where  "c"  is  speed  of  light  and  "h"  is 
Planck's  constant,  it  is  seen  that 


R ft)  = /■»•*  x /d>/9  A Rx  ft) 


(9) 
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Inclusion  of  detector  and  atmospheric  factors  into 
Equations  (2)  and  (3)  results  in  the  following: 


HEFFv  - ^ oU 

heff0  ^ Da)  A(i)  V/yCT^d'X 


(10) 


NEIw 

NEIC 


^■Wv,)D(a)Qx(Vw)^  + fi,Vfr,rw)D(>)tl&<r,)<n 

yd  /t  *\  7^ 


QT)a) 


(11) 


These  are  the  integral  forms  of  the  equations  for  which  a 
numerical  solution  was  developed  and  adapted  to  the  CDC  6600 
Computer.  A fortran  program  listing  appears  in  the  appendix.  The 
solution  in  general  requires  inputs  for  .25  micron  intervals  of 
values  of  window  transmittance  and  emittance,  detector  response, 
atmospheric  radiance  and  transmission.  Values  of  the  blackbody 
function  and  the  photon  quantum  exitance  function  are  calculated 
internally. 
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V - INPUTS 


A large  number  of  variables  is  required  to  be  input  to  LOWTRAN 
U for  any  calculation  of  atmospheric  radiance  and  transmittance 
values.  For  purposes  of  this  paper,  four  model  atmospheres  were 
chosen  to  allow  simple  presentation  of  the  effect  of  different 
atmospheres  on  the  SRF  due  to  a window.  The  term  "atmosphere"  as 
used  here  includes  operational  scenario  factors;  i.e.,  altitude 
and  slant  range,  as  well  as  the  basic  intensive  atmospheric 
properties.  The  four  model  atmospheres  are: 

MOD  1 • ATM  (1)  - 0.30KM  Altitude,  6.06KM  Slant  Range, 

Visibility  3.00KM 

MOD  2 • ATM  (2)  - 6.00KM  Altitude,  9.09KM  Slant  Range, 

Visibility  3.00KM 

MOD  3 • ATM  (3)  - 0.30KM  Altitude,  15.25KM  Slant  Range, 

Visibility  10.00KM 

MOD  U • ATM  (U)  - 6.00KM  Altitude,  15.25KM  Slant  Range, 

Visibility  10.00KM 

All  model  atmospheres  assume  the  midlatitude  winter  model,  a final 
altitude  of  0.03KM  and  a background  temperature  of  300°K.  Trans- 
mittance and  radiance  values  for  these  four  models  are  presented 
in  Table  1. 


Three  different  window  materials  were  considered.  Transmittance 
and  emittance  values  for  each  material  are  presented  in  Figure  2 and 
in  Table  2 and  are  based  on  measured  values  reported  recently  under 
several  government  contracts.  Some  extrapolation  was  necessary  to 
arrive  at  values  for  a window  of  appropriate  thickness  for  PAVE 
TACK.  The  values  are  for  anti-reflection  coated  windows,  where 
reflection  is  either  known  to  be  or  assumed  to  be  on  the  order  of 
1.0  to  2.0  percent  per  surface.  The  values  are  not  considered  to 
be  dependent  on  temperature  (which  is  basically  true);  however, 
the  temperature  of  the  window  is  required  as  an  input  to  the  program 
since  it  is  required  in  calculating  photon  noise  from  the  window. 

Values  for  Detector  Factor,  D(X),  were  obtained  from  Figure  1, 
the  relative  spectral  response  of  the  AN/AAQ-9  FLIR  System,  as  re- 
corded during  laboratory  testing.  Values  of  D(l)  read  from  Figure  1 
in  .25  micron  increments  from  7.0  to  13.0  microns  are  listed  in  Table  2. 


8 


ENA  79-6 


.3 


9**0  10.0 

WAVELENGTH  (MICRONS) 


FIGURE  1 

RELATIVE  RESPONSE  OF  AAQ-9  FLIR 
(Normalized  at  10  Microns) 
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ATMOSPHERIC  TRANSMITTANCE/RADIANCE  VALUES  FOR 
FOUR  MODEL  ATMOSPHERES  IN  .25  MICRON  INTERVALS 
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TRANSMITTANCE  AND  EMITTANCE  OF  THREE  WINDOW  MATERIALS 
AND  AAQ-9  DETECTOR  FACTOR  IN  .25  MICRON  INTERVALS 
MAT  1 - ZnS,  MAT  2 = GaAs,  MAT  3 = ZnSe 
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VI  - RESULTS 


Table  3 presents  the  sensitivity  reduction  factors  calculated 
for  each  of  the  window  materials  (at  300°K)  and  each  of  the  model 
atmospheres.  In  all  cases,  the  SRF  for  zinc  sulfide  is  signifi- 
cantly less  than  for  either  zinc  selenide  or  gallium  arsenide. 

Due  to  the  almost  perfect  transmittance  of  zinc  selenide,  the  SRF 
for  all  conditions  is  almost  constant.  It  should  be  understood, 
however,  that  system  performance  in  terms  of  NET,  MRT  or  other 
performance  measure  is  not  constant  for  different  atmospheric 
models.  A constant  SRF  means  only  that  the  window  degrades  the 
system  performance  by  the  same  factor  regardless  of  the  atmosphere. 

The  values  of  SRF  for  atmospheres  1 and  3,  and  likewise  2 and 
k,  are  nearly  identical  even  for  zinc  sulfide  because  the  spectral 
distribution  of  atmospheric  transmittance  and  radiance  is  very 
similar  for  those  models.  Since  atmospheric  factors  are  present 
with  and  without  the  window,  it  is  the  spectral  distribution,  not 
the  magnitude  of  these  factors,  which  affects  the  calculation  of 
SRF. 

The  effect  of  window  temperature  on  the  SRF  for  different 
window  materials  is  clearly  illustrated  by  Figure  3.  As  tempera- 
ture increases,  SRF  decreases  due  to  increased  noise  coming  from 
the  window.  This  effect  is  independent  of  the  atmospheric  model 
chosen. 


Zinc  Sulfide 

.771 

.698 

.733 

.701 

.733 

Gallium  Arsenide 

.930 

.912 

.922 

.913 

.923 

Zinc  Selenide 

.987 

.986 

.986 

.986 

.987 

TABLE  3 

CALCULATED  SENSITIVITY  REDUCTION  FACTORS 
FOR  WINDOW  MATERIALS  AT  300°K 
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IX  - APPENDIX 


bbst  quaHT*  * 


program  srfact 


74/74  OPT ■ 1 


0P/14/T*  12. 07. Si 


PROGRAM  SRFACT  1 I NPuT  .OUTPUT , TAPESMNPU7. 7AP£6"0UTPC'7I  srf actor 

OIMtNSIOM  QTMI24)  .QTW(P4> .WI24)  SRF  ACTOR 

DIMENSION  0|?O  .A  (?»,  t ) ,M  (?«.«)  SHFACTOR 

— DIMENSION  SHF  13.41  >7(24.])  .£(24.31 SRF  ACTOR 

INTEGER  ATM , ATMAX  SHFACTOR 

RE  »L  Nt  I0.NEIM.NEIR SRFACTOH 

NAMELIST/TRANS/T  SHFACTOR 

NAHEUST/EmISS/E  SHFACTOR 

NAMElIST/OETEC/D  SRFaCTOH 

NAmEliST/ATMOS/A.R  SRFACTOH 

MEW. « * SRFaCTOH 

ML  INC".  PS  . SRFaCTOH 

1 MA  X "PA  SRFACTOH 

T M*  TOO . . ...  SRFaCTOH 

TmOOm.300.  SRFaCTOH 

0EL1 «?5.  SRFaCTOH 

C Till  AI.O  till  ARC  THE  TRANSMISSIVITY  anO  EmISSIVITT  OF  the  MINDOM  FOR  SHFACTOR 
C MAVElCHlilM  INTERVAL  1 1 I THESE  VALUES  MOST  BE  INPUT  SHFACTOR 

C Uldll)  ANO  UlMtl)  ARE  THt  PHOTON  EmITTAnCE  FUNCTION  EVALUATED  FOR  TmC  SRFaCTOH 
C WAVELENGTH  INTERVAL  III  ANO  A I THE  TEMPERATURE  OF  ThE  BACKGROUND  I TBI  OR  SRFaCTOH 
CTEmpehaTURE  OF  The  WINDOW  ITWOOwt.  THE  VALUES  ARE  CALCULATED  BV  ThE  SRFaCTOH 
C SUBROUTINE  PHOTON.  SRFaCTOH 

c win  is  The  GLACKboov  Function  evaluateo  For  The  mavClEngth  interval  ii  srfactor 

C_TEMP£HaIuH£  of  The  BACKGROUND.  SRFaCTOH 

C TB-TEmP  of  bACKGROUNO  SRFaCTOH 

C ml  INC"  wIOTh  OF  each  WAVELENGTH  INCREMENT  IN  MICRONS. . SHFaCTOR 

RLaOIS.Thansi  SRFaCTOH 

HEaOIS.EHISSI  SHFaCTOR 

REaOis.OE IEC)  SHFaCTOR 

REA0(SiATh0SI._ SRFaCTOH 

SRFaCTOH 

ATmaX"!  SRFACTOH 

W M A X" | SHFaCTOR 

IF  I T ( IP, 21  ,NE.0.0)KMAX"2  SRFACTOH 

IF  I r IIP. 31  . NL .O.OIKMAX"]  SRFaCTOH 

IF  I A 1 12.2)  .NE  .0.01  A|MAI"2 SHFaCTOR 

If  (A  ( 1 2, 31  ,NF.  .0.01  AT  MAX"  J SHFACTOR 

IFlA(12.A|.NE.0.0IATMAA>A  SHFaCTOR 

SRFaCTOH 

• .PRINT  INPUT  VALUES  OF  TRANSMITTANCE'  EMITTANCE'  ANO  DETECTOR  FACTOR  SRFaCTOH 

SRFaCTOH 

PRINT  106  SHFaCTOR 

PRINT  loT  SRFaCTOH 

PRINT  I OB  - SRFACTOH 

WLH 1 NaWL 1 SHFaCTOR 

. 00  to  I " I • IMAX  . _ SHFaCTOR 

WLMAX"WlHIN'WL|NC  SRFaCTOH 

PRINT  nO'WLMINtWLMAA'(T(l'KI'Kalt]i<<eillJttJ"i'3l'0<Ii SHFAFIOR 

10  WLM|N"WLMAX  SRFaCTOH 

SHFACTOR 

• PRINT  INPUT  values  of  ATMOSPHERIC  TRANSMITTANCE  AND  RADIANCE  SRFaCTOH 

SRFaCTOH 

PRINT  M2  SRFACTOH 

PRINT  1 14  SHFACTOR 

HLMINaWLI  SHFaCTOR 

00  12  I ■ I • I HAA _ , SHFACTOR 

WLNAXwWLHlN.NLINC  SRFaCTOH 
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MIS  Pi«  IS  BEST  «UU.in  "*»«*** 

S£!iSn«w isa»»»?  — 


PH  I NT  1 16.WLMIN.WLMAA. ( A ( I , J| , J. ( , 4 1 , (R ( l , K | ,K. | , 4) 
12  WLHIN.WLMAX 

DO  23  lltMP.1,5  . 


WLMlN.WLI 

DO  ?S  I*i.|max 

WLMAX*»LHIN.MLINC 


c calculate  values  or  otb  for  each  mlinc.  . if  radiance  values  from 

c LOwTRan  are  input.  ThEST  value's  OF  0TB  WILL  NOT  BE  USED. 

C »LL  pnotoniwlminiwlmaa.u, jmaa.tbi  

OTHIIiaQ 

C CALCULATE  VAIUES  OF  OTW  FOR  EACH  WLINC 

CALL  PHOT  ON  tWLHlN.WLHAX.Q.  JMAX  . T WOQWI 

OTW I 1 1 *0 

c . _ 

C CALCUL  ATI  VAIUES  OF  W FOR  EACH  ML  INC 

CALL  bHODY I WLM I N a WLMAX , wbB , UMAX • TBI 

will . wNB 

2S  WLHINaWLHAX  . 


00  it  AIMM.ATMAX 

IF  IP  I 12, ATM)  . to. 0.0)  l>0  10.29  

(■*00  t « I . I , I HA  X 

^-(WLM  I M>  WL  l)  mHA  * W.  ♦ Wl/EL 

OTRIIIAHII.ATMI  WLMIN.NL MAX  I *.S*1 .SBC19*  .Z&~ 

?B  WLMINaLLMAX  

r.O  TO  30 

29  PRINT*. "VALUES  OF  QTB  ARE  NOT  FROM  LOWTIUN* 


• COMPUTE  HEFFW.  HEFFO.  NEIW.  AND  NCIO 
JO  DO  31  KnI.KMAX 

HEFFW*  0.0  „ .... 

HEFF  0*0 • 0 

NE I w*0 » 0 . 

NE 10*0.0 

00  3S  1*1.1  MA X . . 

HEFFWaH|FFw.TII.KI*AiI,ATM)*DIII*WIII 
hEFF0*n‘FFO.A<I.AT.N)*0(II*Wlll  - _ 

NE I w*NE I w*E II.A)*0ill*0TwlII«TlI.KI*0lll«0TBIII 

JS  NCIO*NEIO. 0(11*018(11  

NE  TR* (NElw/NE 101 **.5 

31  SHF  I K , A T M I w (HEFF W/hEF FOI /N( IM 

27  CONTINUE 


*0  PRINT  140 • TB  t TwOOW  . _ . 

PRINT  US.  ( (K.ATM.SRF  IK.  ATM  I .Kw|  .KNAXl  «ATN*1 , ATMAX) 


PRINT*.  « • 

23  TWDOW.TWOOW.OELT 
to  STOP 


..  |0t  FORMAT 
107  FORMAT 
lot  FORMAT 
CAT  3* . 

110  format 

112  FORMAT 
...  114  FOPmaT 
CX..MOD 
C*./l 

lit  FORMAT 
_ 140  FORMAT 
CN*  I 

_ J4S  - FORMAT 
C 

_ . ..  ENO  _ 


(IMI.40X. ‘TABLE  OF  INPUT  VALUES*!  ■ - 

I//.4RX, ‘TRANSMIT  I ANCF*.2SA,*ENITTANCf*. I J*.*OE T£CTOR*l 
I 9X, .WAVELENGTH  band  (MICMOnsi *.JX.*MAT  I * ,7X. *MAT  2*.7X,.M 
7 X . «MA T I*. IX, .mat  2*,7X.*mAT  3*i7X.* FACTOR* I 
I13X.F7.3.*  TO*.F7.3,7FI2.3l 
I ////.  A 0*  . ‘ATMOSPHERE  FACTORS*,/!' 

( 3X . * wavEL ENOTh  BANO  (MICRONS! *, ‘TRANS  MOO  1*,7X..M00  2*7 
3* , 7 X • *MOD  4***  RAO  MOD  1*,TX.*M00  2.,TX,*M00  J*,7X.*MOO  4 

I1JX.F7.3.*  T0**F7.3,0F12.3! 

(//•10X..TB...FS.1.1 OX. .WINDOW  TEMP.* .Ft. 1 ,21. .DEGREES  xCLVl 


I20X,*SRF  (WINDOW*. 12,*, ATM.  MODEL*. 


SRF ACTOR 

59 

SRFaCTOR 

60  - 

SRFACTOR 

61 

SMFACI0R 

62  — 

SRF  ACIOH 

63 

SRFaCTOR 

66 

SRFaCTOR 

65 

SRFAC10R 

66  ... 

SRF  ACT  OR 

67 

SRFACTOR 

66  . 

SRrACTOH 

69 

SHF  ACTOR 

.10 

SRF  AC  1 OR 

71 

SHFaCIOR 

72  - 

SRFACTOR 

7* 

SRFACTOH 

- 74 

SHFACIOR 

75 

SRFACTOR 

76  . 

SRFACTOR 

77 

SHF  ACTOR 

76  _ 

snr actor 

75 

SRFACTOR 

60  — 

SRFaCTOR 

61 

Srf actor 

. h i _ 

SRFACTOR 

63 

SRFACTOR 

H4  . 

SHF  ACTOR 

6S 

SHF' AC  TOR 

e* 

SIFACTuR 

hi 

SRFACTOH 

he  

SRFACTOR 

h 5 

SRFACTOH 

7* 

SRFACTOR 

^ 1 

SHFACIOR 

92 

SHrACTOM 

73 

SRFaCTOR 

. ... 

SRFaCTOR 

75 

SPF ACTOR 

, * 06 

SRFACTOR 

57 

SRFACTOR 

76 

SRFACTOR 

79 

SRFACTOR 

100  . 

SRF AC  100 

101 

SHFACIOR 

102 

SRFACTOH 

103 

SRFACTOR 

104 

SRFACTOH 

105 

SRFACTOR 

106 

SHF  ACTOR 

10? 

SRFACTOR 

106 

SRFACTOR 

105 

SRFACTOR 

no 

SRFACTOH 

111 

SRFACTOR 

112  _ 

SRFACTOH 

113 

SRFACTOR 

114  „ 

SRFACTOH 

115 

SRFaCTOR 

116 

SRFaCTOR 

11? 

SRFaCTOR 

n» 

SRFACTOH 

119  - 

SRFaCTOR 

120 

SRFaCTOR 

121  - 

SRFaCTOR 

122 

SRI  ACTOR 

123  - 

SRFACTOH 

124 

SRFACTOH 

125  * 

SRFaCTOR 

126 

SRFaCTOR 

12?-- 

SRFaCTOR 

120 

SRI  AC I OR 

129 

SRFaCTOR 

130 

SRFaCTOR 

131  - 

SRFaCTOR 

132 

SRFACTOH 

-til- 

SRFaCTOR 

134 

SRFaCTOR 

115  . 

P 


® jo  sec 


SUBROUTINE  PHOTON 


T4/T4  OPT  * 1 


It/U/ll  IT.tl.ll 


SHI  AFT  OH 
SHF AC10M 
SHE  PC TOM 
sofactor 
Srfacior 
srf  »c‘oh 
Su»  »cm» 
SRFacTom 
swr actom 
SUFtCTOfc 
ShFAL'Ioh 
SHF  ACTO* 
SHFACTOH 
SHFACTOH 
SHFaCToh 
SBFaCIoH 
SHF  AC  I OH 
SHF ACTON 
SxFaCTuH 
SHF  ACT  OH 
SHF' ACT  OH 

SHFaCToh 
SHFACTOH 
SHFACTOH 
SHF  ACT  Oh 


subroutine  photonihlmin.hlmax.o.jmaa. TEMPI 

ALL  OMITS  of  LENGTH  AHE  M1CHONS,  TEMP  is  KELVIN 

OELHLA.OI  - 

PIO.I4I6 

0*2.997  76E14  ...  . 

C2*l .43BE4 

UMAX* IWLHAX-WLM1NI/DCLML..S 

ML«HLMIN  , — 

C*?.*PI*CI 

QQlO«C/  |*L**4.<  IEXPIC2/(*LMEHP|  IM.IL 

00  20  JM.JNAX 

hl-al-oelhl  

QnF**C/1*LA*4.a(EXPIC2/I*LMEMP»»M,I» 

0*0.  ,s»  <une*.qoloi  •uclml  _ 

ogloaqnem 

I CONTINUE  — - 

THE  units  OF  a ARE  NO*  CHANGED  FROM  PHOTONS  • NlCR0N**-2»SEC**-l  TO 
_.  PHOTONS  • CM»«-2  • SEC»»-l  


TA/lA  opt«i 


02/1  A/T*  I2.07.SS 


SUBROUTINE  BBOOT IHLMIN.ULNAX.MBBf JMAXf TEMPI  SHFACTOH  161 

ALL  UNITS  OF  LENGTH  ARE  MICRONS.  TEMP  IS  KELVIN SRFACTOR  162 

OELHLA.OI  SHFACTOH  |A3 

PI-3.1416  _ SMFACTOR  ISA 

Cl *2. 99  7 76CI A SHFACTOH  US 

. . C2M.43F.EA  _ SHFACTOH  166 

PM«6.6?bC-t!  SHFACTOH  1 AT 

C*2.«PI»C1A.*.*PH  ....  SRFACTOR  lbB 

*88*0.0  SMFACTOR  I A A 

JMAXA IHLMAX-WLMINI/OELVL..S  , SRFACTOR  1 To 

NL-WLMIN  SHFACTOH  1>| 

_ *OLO*C/  l*L**S.A  (CAP  IC2/  IHLMENPI 1 Ms  I 1 SHFaCI0h_JT2 

00  20  JM.JMAK  SHFACTOH  1T3 

*L*  *L . OELHL  SHFACTOH  | TA 

*N(.*C/I*L*aS.*IEXPIC2/IHL*TEHP))-|,)I  SHFACTOH  ITS 

--  *8«*«BB..Sa|BNCR**0l01»0EL»*L SHF  ACT  OH  | TA 

*OLO*hNEV  SHFaCToh  ] TT 

0 CONTINUE  SHFaCToh  ) TB 

The  units  of  MSB  ARC  NO*  CHANGED  FROM  ERG/ (SEC*M|CR0Na*2I  TO  MATTS/CM  SHFaCTOR  |!« 

. . .... SHFaCToh  100 

VBBaUBBMB.  SHFACTOH  181 

SHFACTOH  182 

RETURN  SMFACTOR  IB] 

CMO — SRFACTOR  -4  BA 
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